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ABSTRACT 

Since launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has de¬ 
tected many hundreds of bursts from magnetar sources. While the vast majority 
of these bursts have been attributed to several known magnetars, there is also 
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a small sample of magnetar-like bursts of unknown origin. Here we present the 
Fermz/GBM magnetar catalog, giving the results of the temporal and spectral 
analyses of 440 magnetar bursts with high temporal and spectral resolution. 

This catalog covers the first five years of GBM magnetar observations, from 
July 2008 to June 2013. We provide durations, spectral parameters for various 
models, fluences and peak fluxes for all the bursts, as well as a detailed tem¬ 
poral analysis for SGR J1550—5418 bursts. Finally, we suggest that some of 
the bursts of unknown origin are associated with the newly discovered magnetar 
3XMM J185246.6 + 0033.7. 

Subject headings: catalogs — pulsars: individual (SGR J1550—5418, SGR J0501-I-4516, 
1E1841-045 , SGRJ0418+5729, SGR1806-20, SGRJ1822.3-1606, AXP4U 
0142+61, AXPIE 2259+586, 3XMMJ185246.6 + 0033.7) - stars: neutron - 
X-rays: bursts 


Introduction 


In 1979 Mazets et ah (1979a|b| reported the detection with the Venera satellites of a 


series of short, energetic events apparently originating from the same directions on the sky. 


Although their durations would place them in the short gamma-ray burst class, Mazets et 


al. (1979a) argued that the recurrent nature of the emission indicated that we were looking 


at a new phenomenon. In the following decade, these sources were indeed recognized as 
a different class of objects, called soft gamma repeaters (SGRs), based on their repeated 


high-energy burst activity (Golenetskii et al. 1984, Laros et al. 1987; Atteia et al. 1987 


Kouveliotou et al^|1987 ). In the late 1990s it was discovered that these objects were isolated 
neutron stars characterized by extremely strong dipole magnetic helds of 1014 _ iq15 g 
(K ouveliotou et al.|1998 1999), right along the lines of their theoretical prediction by Duncan 


& Thompson (1992). Gurrently, the bulk of the magnetar population comprises mainly two 


classes of objects: SGRs and Anomalous X-Ray Pulsars (AXPs). Their distinguishing traits 
are slow spin periods (~ 2 — 12 s), large period derivatives (~ 10“^^ — 10“^° s s“^), and 
periods of activity during which they emit multiple bursts in the hard X-ray/soft y-ray 
energy range. With respect to the latter burst characteristic, magnetars seem to fall into 
two groups: sources that emit a few hundreds to a thousand bursts per active episode and 
sources emitting only a handful of bursts when active (e.g., Gogu§j20'l4 ). For detailed: 


. reviews 


on magnetars, see Woods & Thompson (2006) and Mereghetti (2008) 


The magnetar discovery rate increased signihcantly when wide £eld-of-view monitors 
with location capabilities became operational, or multiple instruments in space allowed for 
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triangulation of their signals to be performed. The hrst wave of discoveries in the late 80s 
was enabled mostly by the Vela and Prognoz satellites. This was followed by the synergistic 
era of the Burst And Transient Source Experiment (BATSE) onboard the Compton Gamma 
Ray Observatory, the Rossi X-ray Timing Explorer (RXTE), and BeppoSaX in the 90s; and 
hnally, in the mid to late 2000s, the start of the era of NASA’s Swift and Eermi observatories. 
The latter have onboard detectors dedicated to monitoring the high-energy transient sky 
(the Swift/Burst Alert Telescope (BAT) and the Fermi/Gamma-ray Burst Monitor (GBM), 
respectively), and have been extremely efficient in complementing each other in detecting 
and identifying new magnetar sources. 


Since its activation in July 2008, GBM has recorded several hundreds of bursts from new 
and known magnetars, and a few magnetar-like bursts from unknown sources within rather 
large error box regions. A large number of papers has already been published by the GBM 


magnetar team on some of these sources with multiple bursts (SGR J1550—5418: 

Kaneko et 

al. 

2010 

Lin et al. 

2012 

van der Horst et al. 2012 von Kienlin et al.||2012l Huppenkothen 


et al.||20T4 ; Younes et ah 2014 SGR J0501-I-4516: Watts et ah 2010 Lin et ah] 2011a 


Huppenkothen et al. 2013; SGR J0418-I-5729: van der Horst et ah 2010). Here we present a 


comprehensive overview for those magnetars, including bursts that have not been analyzed 
before, and also for magnetars with only a few events (SGR 1806—20, SGR J1822.3—1606, 
AXP4U 0142-1-61, and AXP IE 2259-1-586). We have compiled the entire GBM magnetar 
burst sample in a catalog for the hrst hve years of the mission. In Section we describe 
the GBM data types, selection criteria, and analysis methods. We discuss the properties of 
bursts from conhrmed magnetar sources in Section]^ In Section]^ we discuss bursts for which 
we do not have conhrmed source identihcations. We conclude in Section with a comparison 
across the source properties and a discussion on the characteristics of all magnetar bursts 
detected with GBM. 


2. Data Analysis 


GBM is an all-sky monitor with an 8 sr held of view, consisting of 12 Nal detectors, 
covering an energy range from 8 keV to 1 MeV, and 2 BGO detectors, covering 0.2 — 40 MeV 


(Meegan et al. 2009). GBM has three data types, of which two (GTIME and GSPEG) are 
continuously recorded and one (time-tagged event data, or TTE) recorded only after GBM 
had triggered on a burst or other selected time intervals during most of the catalog period. 
This conhguration changed after November 2012, and currently GBM is recording TTE 
data continuously. While GSPEG data have good spectral resolution (128 energy channels 
over the full energy range), its temporal resolution of 1.024 s is too long for the ~ 0.1 s 
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long magnetar bursts. The temporal resolution of CTIME data (64 ms) is better suited 
in this respect, but the spectral resolution (8 energy channels) is not sufficient to perform 
detailed spectral analysis. TTE data have the highest temporal (2 /is) and energy (128 
energy channels) resolution, and are therefore perfectly suited for magnetar burst analyses. 
The downside is that this data type is typically only available from ~ 30 s before a burst 
trigger to ~ 300 s afterwards, followed by another ~ 300 during which no new trigger can 
take place. For further details on GBM and its data products, see Meegan et ah (2009). 


2.1. Data Selection 

Given the spectral softness of magnetar bursts, with typically no emission above 200 keV, 
we only use Nal data for our analyses. The 12 Nal detectors are oriented in such a way that 
GBM can monitor the entire sky, which implies that only a subset of the detectors can 
see one specihc location on the sky at any given time. For every burst we have used only 
detectors with an angle between the source direction and the detector normal of less than 
60°, to ensure a large effective area for these events. We note that in previous studies of 
GBM bursts the viewing angle varies from 40 to 60°, but this does not have a signihcant 
effect on the temporal and spectral results. We have excluded detectors when a source was 
obstructed by parts of the spacecraft or the Large Area Telescope (LAT) onboard Fermi. 

We have identified a total of 427 triggered bursts associated with magnetars between 
July 2008 and June 2013, of which the associations have been established based on local¬ 
ization, detections by other instruments, or source bursting activity. We have also searched 
for, and included, all untriggered bursts that happened during the 30 s before and 300 s 
after each trigger, since TTE data were available at these times. Finally, we analyzed 19 
magnetar-like bursts during this period, i.e. bursts with durations and spectra similar to 
those observed in magnetar bursts. Due to the poor positional accuracy of GBM (in par¬ 
ticular for spectrally soft events), the crowded areas on the sky where they were found, and 
the lack of triggers by other gamma-ray instruments, it was not possible to confirm their 
nature. Several other untriggered bursts were found during the catalog time period, but are 
not included here since they lack TTE data. In Table [T] we show the magnetars included in 
this catalog, their burst active periods, and the total number of bursts per magnetar used 
in our analyses. 
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2.2. Temporal &: Spectral Analysis 


The majority of bursts presented in this catalog have been published (van der Horst et 


al.|2010 ; Lin et all2011a|[b 2012 van der Horst et ahj2012 ; von Kienlin et al.||2012 Younes et 
al. 2014). For the remaining bursts we performed temporal and spectral analyses following 
the techniques and criteria as described in detail in these publications. The Tgo and T 50 
durations, dehned by the time that the cumulative counts rise from 5% to 95% and 25% to 


75%, respectively (Kouveliotou et al. 1993), were calculated in photon space, i.e. by using 


the intrinsic (deconvolved) burst spectra instead of the detector recorded counts to dehne 


the durations (for a detailed description, see Lin et al. 2011a). Since many of the bursts 
consist of multiple peaks, we required for two events to qualify as two separate bursts if the 
time between their peaks is longer than a quarter of the magnetar spin period, and the count 
rate drops to the background level between the peaks. We excluded the brightest bursts that 
saturated the High Speed Science Data Bus of GBM from the duration calculations, and for 
these bursts we only used the non-saturated parts in our spectral analysis. The spectral 
analyses were performed with the software package {v4-3), using detector response 

matrices generated with GBMRSP {v2.0). To generate the response matrices we used the best 
available positions from other satellites or telescopes for the magnetar bursts in Section 
and the GBM positions for the bursts of unconhrmed origin in Section]^ We htted the time- 
integrated spectra over an 8 — 200 keV energy range to different spectral models: a single 
power law (PL), a blackbody (BB), a Gomptonized model (GOMP), an optically thin thermal 
bremsstrahlung model (OTTB), and two blackbody functions (BB-I-BB). In each case we 
used the Gastor G-statistic to determine the model signihcance. For SGR J1550—5418 we also 
performed time-resolved spectroscopy of the bursts at 4 ms time resolution to determine their 
peak fluxes. For this purpose we only used the GOMP model. Finally, temporal analyses 
were performed on the whole sample of SGR J1550—5418 bursts, and the methodology and 
results are discussed in Sections 13.1.41 and 13.1.51 


3. Magnetar bursts from confirmed sources 


3.1. SGRJ1550-5418 


SGR J1550-5418 (RA = 15^50“54M1, Dec = -54°18'23".7; Gamilo et al. 


became 


active in October 2008, when it emitted several hundreds of magnetar-like bursts (Palmer 


^R.S. Mallozzi, R.D. Preece, & M.S. Briggs, ’’RMFIT, A Lightcurve and Spectral Analysis Tool ©2008 
Robert D. Preece, University of Alabama in Huntsville, 2008 
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2009 Kouveliotou et al. 2009); this activity ceased April 2009. The source was discovered 


by the Einstein satellite (Lamb & Markert 1981), named IE 1547.0—5408, and identihed as 


a magnetar candidate following observations with XMM-Newton and the Chandra X-Ray 


Observatory (Gelfand & Gaensler 2007). It was later confirmed as a magnetar with a spin 
period of 2.07 s, a period derivative of 2.32 x 10“^^ s/s, and an inferred magnetic field of 


2.2 X 10^“^ G (Gamilo et al. 2007 for a detailed description of the history of SGR J1550—5418, 
van der Horst et al. 2012 and von Kienlin et al. 2012). Although SGRJ1550—5418 


see 


was initially classihed as an AXP, after this extremely prolihc activity very similar to ones 
from SGRs 1806—20, 1900-1-14, and 1627—41, it was reclassihed as an SGR and renamed 


SGRJ1550-5418 ( 

Palmer 

2009; 

Kouveliotou et al. 

2009D. The 

source distance has been 

estimated to be ~ 

5 kpc (( 

Jelfand & Gaensler 

2007 

Tiengo et al. 

2010 

)• 


Throughout this section, we will refer to the activity of SGR J1550—5418 in terms of 


three active periods: October 2008 (von Kienlin et al. 2012), January 2009 (van der Horst 


et al. 

2012 

), and March - April 2009 ( 

von Kienlin et al. 

2012 

)• 

von Kienlin et al. 

(2012) 


studied the GBM burst data for the first and third periods and found that the spectra of 
the former were best fit with a single BB, while the latter were best fit with an OTTB. 


van der Horst et al. (2012) studied the bursts of the hrst two days of the second period. 


and found that their spectra were £t equally well with a GOMP, an OTTB, and a BBj-BB 
model. Lin et al. (2012) later refined the results of van der Horst et al. ( 2012| ), using bursts 
from the second period that were also seen with the Swift/X.-'Ray Telescope (XRT), and 
established that a BBj-BB model best £t these joint spectra. Taken together these three 
studies suggest that the first outburst was significantly different than the other two periods. 


which exhibited similar spectral and temporal properties. Younes et al. (2014) presented 


time-resolved spectroscopy for the brightest bursts over all three periods, and obtained an 
estimate of the lower limit on the internal magnetic field. Besides the analysis of these bursts. 


Kaneko et al. (2010) presented spectral and temporal analysis of enhanced persistent emission 


during the onset of the second bursting period, resulting in the discovery of the smallest hot 
spot ever measured for a magnetar. Furthermore, recent detailed variability analysis has 
revealed candidate quasi-periodic oscillations (QPOs) in bursts during the second emission 


period (Huppenkothen et al. 2014) 


We searched for additional bursts in the GBM data during the gaps between the pre¬ 
viously studied bursting periods, and found 66 previously unstudied bursts that temporally 
best belong to the second period (spanning 22 January — 24 February 2009). These 66 
bursts include 5 untriggered bursts for which we have TTE data, which were found with 


the search algorithm described in Kaneko et al. (2010). Between 24 February and 22 March 


there were no bursts observed from SGR J1550—5418, although the source was still visible 


in the X-rays (see e.g., Scholz & Kaspi 2011) 
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3.1.1. Temporal and Spectral Evolution 


We study here the statistical properties of all bursts from SGR J1550—5418 detected 
with GBM and the evolution of the burst activity across the entire 2008 — 2009 source 
activation period. Figure [^displays the Tgo and T50 duration distributions of 354 unsaturated 
events (see also Table [^. We fit these distributions with a log-normal function (solid line) 
and find that the ^90 (T50) duration centers at 155 ± 10 ms (51 ± 5 ms), with a logarithmic 
width of ~ 0.4. The duration distribution of the bursts from SGR J1550—5418 is similar to 
what we have seen in other magnetars (e.g., Woods et al.||1999 


et ah 2001 Gavriil 


et ah 2004 Esposito et ah 2008; Lin et ah 2011a[ ). Figure shows the evolution of Tgo and 
T50 over time. We hnd no significant trend in the evolution of the burst durations over the 
entire outburst period of roughly 7 months. We note that the longest event on Figure is 
very faint with a peak flux of 9.2 x 10“^ erg/cm^/s (see also Lin et ah 2013 for a sample of 


very dim but long SGR J1550—5418 bursts detected by Swift/XKT). 

In Figures |3]-|^ we present the evolution of time-integrated spectral parameters from 
the OTTB, GOMP and BB-I-BB models for the entire sample of SGR J1550—5418 bursts 


van der Horst et ah 

2012; 

von Kienlin et ah 

2012 


temperatures of our BB-I-BB fits for the second and third periods along with the temperature 
of the single BB that best fit the first activity period. The latter is included in both panels 
(high and low kT) for comparison purposes. As is evident from the lower panel (high kT), the 
temperature of the single BB is similar to the high kT of the BB-I-BB model, which means 
that the lower kT BB is absent in the hrst activity period. This conhrms the suggestion 


of von Kienlin et ah (2012) that there is a spectral evolution from the hrst to the other 


burst active periods. We note that it would be possible for the lower kT BB to be present 
below the detection limit of GBM, but this would still indicate a signihcant change in the 
spectral shape between the different periods. The burst peak energies, derived with the 
GOMP (Tpeak) and OTTB (kT) models, are presented in Figure]^ in the upper and lower 
panel, respectively. We observe no signihcant evolution across the outburst. However, the 
evolution of the spectral index of the GOMP model (Figure conhrms that the bursts of 
the hrst burst period are signihcantly harder. 


The distributions of the spectral parameters discussed above are shown in Figures |6]-[^ 
The solid lines represent Gaussian hts, giving similar results to our prior analyses of indi¬ 
vidual periods (Table |^. We hnd that the BB-I-BB low and high temperatures center at 
~ 4.6 keV and ~ 15.0 keV, respectively, with the width of the higher kT distribution about 
three times broader than the lower kT one. The apparent sub-peak around 12 kT for the 
high-kT distribution is caused largely by bursts from the October activity period. The peak 
energies (Figure]^ from the GOMP and OTTB models agree with each other, with an av- 











































erage of ~ 40 keV. Finally, we find the COMP power-law index distribntion to be narrow 
with an average of —0.93 ± 0.02 (Fignre|^. The signihcant tail excesses in the OTTB kT 
and the COMP index distribntions, above ~ 80 keV and 0, respectively, are dne to events 
from all three periods. 


3.1.2. Energetics 


The flnence range of the SGR J1550—5418 bursts covers over three orders of magnitude, 
from ~7 X 10“® erg/cm^ to at least ~1 x 10“® erg/cm^ (8 — 200 keV). The corresponding 
total energy range is ~ 2 x 10^’^ erg to ~ 3 x 10^° dl erg (with d^ the distance to the 
source divided by its estimated distance of 5 kpc). Compared to the other prolific magnetars, 
these bursts are brighter than the ones reported for SGR 1806—20 (1.2 x 10“^*^ — 1.9 x 10“^ 


erg/cm^; 3.0 x 10^® — 4.9 x 10^® erg; for E > 25 keV; Gogu§ et ah 2001) and SGR 1900-1-14 
(1.2 X 10“^® — 3.3 X 10“’^ erg/cm^; 7 x 10^® — 2 x 10^® erg; for E > 25 keV; Gogu§ et ah 


2001). Woods et ah (1999) give the peak luminosity range of SGR 1627—41 bursts to be 


1039 _ erg/s {E > 25 keV), which corresponds (assuming a burst average duration of 
0.1 s) to an energy range of 10^® —10^^ erg, indicating that these bursts were at the higher end 
of the energy distribution. However, since these data were obtained with a higher threshold 
instrument with limited trigger algorithm options for very short events (CGRO/BATSE), the 
non-detection of a fainter subset may well have been an instrumental effect. SGR J1550—5418 
bursts are an order of magnitude fainter than those of SGR J0501-I-4516 (4 x 10“® — 2 x 10“® 
erg/cm^; 2.0 x 10^^ — 1.0 x 10^® erg; for 8 — 200 keV; Lin et ah 2011a). Finally, all SGR 


bursts are significantly more energetic when compared to bursts from AXP 1E2259 -|- 586, 


which range between 5 x 10^^ and 7 x 10^® erg (2 — 60 keV; Gavriil et ah 2004) 


Next we derive a lower limit on the total energy flnence emitted by bursts from the 
source. The limit is taking into account that we could not include bursts that occurred when 
Eermi was in the South Atlantic Anomaly, the Earth occulted the source, or no TTE data 
existed. For this limit we include the unsaturated parts of the saturated bursts across all 
outburst periods. We get a value of 2.2 x 10“^ erg/cm^, which corresponds to a (lower-limit) 
energy release of 6.6 x 10^^ d\ erg (8 — 200 keV). This is 4 orders of magnitude above the 
upper limit derived for bursts from AXP 1E2259 -|- 586 (2 — 60 keV; Gavriil et ah [2004 ). 

In Figure we present a Log N(>S) — Log (S) diagram for all 354 unsaturated events, 
where S is the energy flnence (8 — 200 keV). We ht this histogram with both a PL and a 
broken PL. For the single PL we only £t the bins >1 x 10“^ erg/cm^, obtaining an index of 
—0.9 ± 0.1. We find that the broken PL fits the data best with a break at 1 x 10“^ erg/cm^, 
and indices of —0.14 ± 0.04 and —0.9 ± 0.1 below and above the break, respectively. The 
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PL slope above the break is steeper than what has been fonnd for other magnetars with a 


large sample of bursts, ranging from —0.5 to —0.7 ( 

Gogu§ et al. 

1999; 

Woods et al. 

1999 

Aptekar et al.||2001 

Gogu§ et al. 

12001 

Gavriil et al.||2004 

), although a similar slope has been 

found for a large sample of SGR 1806—20 bursts ( 

Gotz et al. 

2006 

. The PL slope below the 


break is clearly much shallower, which is due to the lower sensitivity of the instrument to 
faint bursts. Regarding the PL slope at high fluences, the energy and fluence ranges that are 
probed in our SGR J1550—5418 sample differ from those in other magnetar samples, making 
a clear direct comparison unfeasible. 


Finally, we find that the 4 ms peak fluxes of SGR J1550—5418 bursts range from ~8 x 
10“^ to ~2 X 10“^ erg/s/cm^, corresponding to a peak luminosity range of ~5 x 10^® to 
~1 X 10^^ erg/s. The 4 ms fluxes extend the peak flux range to about a full order of 
magnitude higher than was covered with the average fluxes from this source (Lin et ah]|2012 


van der Horst et al. 2012 von Kienlin et ah 2012). Figure 10 shows the evolution of peak 


fluxes across the outburst and demonstrates the uniformly low peak fluxes of the events of 
the first period. 


3.1.3. Correlations 


We searched for correlations among burst properties across the entire burst activity 
from the source, using the Spearman rank order correlation test. We considered multiple 
combinations of parameters and estimated the Spearman coefficient and probability for each; 
the results are listed in Table[^ In most cases, the significance of the correlation was very low, 
but we identified the following significant correlations: fluence — peak flux, i^peak — fluence. 


and F^peak — peak flux (all determined using the GOMP model). Figure 11 exhibits the 
strong correlation between fluence and peak flux; the data are £t with a PL with an index 
1.2 ± 0.4. The F^peak ~ fluence correlation is moderately weaker than the one previously 


reported for bursts of only the second period (van der Horst et al. 2012), but there still is 
a very small chance probability of 9.7 x 10“^^. We previously found a marginal correlation 
between Tgo and fluence (van der Horst et al. 2012); here we find a weaker Tgo — Fluence 
correlation (chance probability of 2.5 x 10“"^), and no correlation between Tgo and peak 


flux (Figure 12). However, we do find a marginal anti-correlation (chance probability of 


7.1 X 10“®) between T 50 and peak flux (see also Figure 13). This can be explained by the 
fact that for bursts of which the light curve is dominated by one very bright peak, the T 50 
is likely to encompass only that bright peak and not include any lower level activity that 
is part of that burst (while the latter is included in the Tgo). This is further illustrated in 
Figure 14, where we show a histogram of the Tgo / T 50 ratio. This histogram peaks between 
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2 and 3, but has a broad tail extending above 10, which shows that there is indeed a large 
amount of bursts that consist of a bright peak and relatively long time span of lower level 
emission. 


We estimated peak fluxes at a 4 ms time scale, and investigated the E-^eak — peak 
flux relation at this time scale for the first time in this study (Figure [l5|. We find that this 
correlation is considerably stronger than the one using the average fluxes, reported previously 
by van der Horst et ah (2012). We £t the data to both a PL and a broken PL. The PL 
index is —0.06 ± 0.01, and the broken PL indices are —0.23 ± 0.03 and 0.16 ± 0.03 with the 
break at 2.7 x 10“^ erg/s/cm^. The break value is almost an order of magnitude higher than 
the one estimated using the average flux values (van der Horst et ah 2012), but closer to 
the break flux determined in time-resolved spectroscopy of the brightest SGR J1550—5418 
bursts (Younes et al. 2014). The latter is also comparable to the break flux estimated for 
the time-resolved spectroscopy of bright bursts from SGR J0501-I-4516 (Lin et al. 2011a). 
We note that we fi; 

(Spearman rank = 


We note that we find a considerably stronger correlation among the bursts below the break 
-0.422, chance probability of 5.9 x 10“^^ 


than the one reported in 


van 


der Horst et al. (2012). For the bursts above the break, we again see an improved but still 


insignificant correlation (chance probability of 0.17) as compared to one found when using 
the average flux (chance probability of 0.64). 


3.1.4- Burst Peak Arrival Times - Pulse Phase Correlation Analysis 


We have studied the burst peak arrival times with respect to the spin phase of the source. 
To determine a burst peak, we fit for each burst the time-resolved spectra accumulated over 
4 ms bins to a GOMP model using RMFIT. The binning process was repeated for different 
bin edges to account for windowing effects. We assigned the time of the center of the brightest 
4 ms as the peak in each burst. All times were barycenter corrected using the JPL ephemeris 
hie as provided by HEASOFT-6.11.1; the phase of each burst peak was estimated using the 


spin ephemeris of Dib et al. (2012). 


Figure 16 shows a histogram of all burst peak phases, with the SGR J1550—5418 pulse 


prohle as measured with XMM-Newton (Lin et al. 2012) overplotted. We hnd that the peak 


arrival times are distributed uniformly across the spin phase (average bin height: N = 17 
events, <7 = 5 events). We tested for a correlation between the peak hux, huence, and burst 
duration against spin phase. Figure [T7| shows the individual points, their weighted averages 
by 0.05 phase bins, and the SGR J1550—5418 pulse prohle. We hnd no correlation between 
these three burst properties with spin phase. 
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3.1.5. Burst Profile Epoch Folding 


Scholz & Kaspi ( 2011[ ) found, using Swift/X.KT data (0.2 — 10 keV), that the epoch folded 


light curves had one (fairly smooth) peak, at 0.45 — 0.65 of the pulse phase. However, this 


was not found in an analysis of Swift/XBUi bursts by Lin et ah (2012). We have performed 


a similar analysis with the GBM data to check this at higher energies. We epoch folded the 


background-subtracted light curves of each burst using the ephemeris of Dib et al. (2012). 


We then added all folded light curves together to create the pulse profiles shown in Figure [T^ 
We display the folded profiles in four energy bands: 8 — 20 keV, 20 — 50 keV, 50 — 200 keV, 
and 8 — 200 keV together with the pulse profile in 0.5 — 10 keV. We notice a peak at phase 
0.15 — 0.20, which appears in all energy bands and is not associated with the pulse maximum 
or minimum. 

We have studied the significance of this peak by calculating its deviation from the mean 
as (Tj = {xmax — x) /o'lg, where Xmax is the value of the peak, and x is the average of the 
histogram without that maximum bin. For the 8 — 200 keV energy range, we find cxj ~ 
5.9. To test for the significance of this deviation, we performed a Monte Carlo simulation, 
randomly shifting all SGR J1550—5418 light curves in phase space by a random number 
chosen from a uniform distribution. After each shift, we created a new count-phase histogram 
and recalculated cTj. This process was repeated 10,000 times, and we find that 1.24% of the 
time, CTj would be at least as high as it is in the real data. This corresponds to a ~ 2.5cr 
result, rendering this feature insignihcant. Interestingly, we hnd that if we remove the top 
ten brightest bursts from the histogram, the peak disappears (Figure [^, suggesting that 
the brightest bursts are preferring a certain spin phase, a fact that may have been buried in 
the arrival time analysis. We tested this peak in the same way as described above and find 
that in this case, cTj ~ 11.4. However, a Monte Carlo simulation on the brightest 10 bursts 
shows that this peak is only a ~ 2.95 ct result. Therefore, our findings confirm the results of 
Lin et ah] (2012) that there is no strong evidence for a correlation between the burst counts 
and spin phase. 


3.2. SGRJ0501+4516 


SGR J0501-I-4516 was discovered with Swift/BAT (Holland et al. 2008 Barthelmy et ah 


2008) when it emitted a magnetar-like burst on 2008 August 22, which also triggered GBM 


(Lin et ah]|2011a). The entire source outburst lasted for about 2 weeks, with multiple bursts 


observed with Swift, GBM, RATE, Konus-Wind, and Suzaku (Enoto et al. 2009 Aptekar et 


aL]|2009 Kumar et al.||2010 ; Lin et al.]|2011a[ ). Gogu§ et al. (2008) found the spin-period with 
RATE: ~ 5.76s; the combined RATEand Swift/X.BT data revealed a spin-down rate of P = 
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1.5 X 10 s/s, which corresponds to a magnetic field of 2.0 x 10^^ G (Woods et ah 2008 


Rea et al. 

2009; 

Gogu§ et al. 

2010 ) 

found with 
0 ".ll (la; 

Chandra to be RA = 05' 

Gogu§ et al.||2010 

). The s 


The most accurate location of SGR J0501+4516 was 


2 kpc (Xu et 


al. 


Lin et al. 2011a). 


GBM triggered on 26 bursts during the entire outburst interval; an untriggered search as 


described in Kaneko et al. (2010) yielded 3 more bursts with TTE data, resulting in a total 


of 29 bursts. Two of these bursts saturated the High Speed Science Data Bus of GBM and 


were excluded from the detailed temporal and spectral analyses which are described in |Lin 
et al. (2011a). In Tables [ 6 ]-[^ we provide the Tgo and T 50 durations in photon space, and the 


spectral fitting parameters for these bursts. We display in Figure 20 the Log N(>S) — Log 


(S) diagram for SGR J0501+4516bursts, computed in this study. We fit the histogram with 
both a PL and a broken PL. For the former we only fit the bins > 1 x 10“^ erg/cm^, yielding 
an index of —0.7 ± 0.1. The broken PL indices are —0.30 ± 0.07 and —0.73 ± 0.08 below and 
above a break at 10“^ erg/cm^, respectively. We note that in contrast to SGR J1550—5418, 
the single PL index and the high-fiuence index for the broken PL are consistent with indices 
found for other magnetars (Section 3.1.2 Gogu§ et al.||1999 ; Woods et al.||1999 


et al. 


2001 Gavriil et al. 2004) 


Finally, we note that a detailed variability analysis of the GBM bursts of SGR J0501+4516 


did not result in a significant detection of QPOs (Huppenkothen et al. 2013) 


3.3. IE 1841-045 


IE 1841—045 is a magnetar located at RA = 18’^41™19®.343, Dec = —04°56'11".16 (from 
Chandra observations Wachter et al.||2004 ), and at a distance of ~ 8.5]'^};o kpc ( |Tian fc Leahy 
2008). ASCA observations have revealed a period of ~ 11.8 s (Vasisht & Gotthelf 1997); 
the spin down rate of P = 4.16 x 10“^^ s/s has been determined through ASCA, RXTE, 


and BeppoSaX observations (Gotthelf et al. 1999, 2002). Swift/BAT was first triggered 


by one burst on 6 May 2010 (Kumar & Safi-Harb 2010). This was followed by a period of 


activity between 8 February and 2 July 2011. During this time. Swift observed 4 bursts (e.g.. 


Barthelmy et al. 2011; Rowlinson et al.| 2011 Melandri et al. 2011). Fermi/GBM triggered 
on magnetar-like events six times during this time interval, with one of these triggers being 
one of the four Swift detections. The other five detections have no concurrent observations 
with other instruments. In addition, an untriggered burst search was performed over all of 


February 2011 and 10 June — 6 July 2011 using the algorithm described in Kaneko et al. 
(2010). This search yielded only one burst, but one for which there is no TTE data, so it 





















































































13 


is excluded from analyses. We argue that despite the lack of concurrent observations, given 
the coincidence in time and position of these bursts, it is highly likely that they originate 
from IE 1841—045 (see 


Lin et al. 2011b 


Detailed spectral and and temporal analyses of these GBM and Swift bursts are de¬ 
scribed in Lin et al. (2011b). Here we provide in Tables [ 8 ]-[^ the Tgo and T 50 durations in 
photon space, and the spectral parameters for these bursts. New in this catalog compared 
to previous studies, we provide OTTB time-integrated spectral hts and 4 ms peak fluxes for 
these bursts (see Table . 


3.4. Bursts from other known magnetar sources 


In this section, we cover 7 bursts each originating from a known source as conhrmed with 
Swift observations, but with not enough bursts to perform a statistical analysis as with the 
magnetars in the previous sections. These include two bursts from a new magnetar source, 
SGR J0418-I-5729, discovered with GBM and also observed with Swift/BAT and Konus-RF 


van der Horst et al. 

2009 

); the GBM data were published in van der Horst et al. 

( 2010 ) 


and are also included in Tables 10-11 The other events originated from SGR 1806—20 


(contemporaneous Swift/BAT detection; 

Bhat et al. 

2010 

); SGR J1822.3—1606, which is a 

new source discovered and localized by Swift ( 

Gummings et aL||2011 

Pagani et al.||2011i 

von 


Kienlin et al. 2011); AXP4U 0142-1-61, for which GBM observed a burst simultaneous with 


a Swift/BAT burst (Palmer 2013 Bhat et al. 2011); AXP IE 2259-1-586, for which a single 
burst was observed with a location consistent with the source during an increased activity of 
the source’s persistent emission (Foley et al. 2012 Archibald et al.| 2013); SGR J1550—5418, 


which emitted a single burst in June 2013, which was also detected with Swift/BAT (Holland 


et al.||2013 ), and was comparable in brightness to the bursts in the October 2008 period. For 
all bursts we performed spectral and temporal analyses following the methods described in 
Section]^ The results of these analyses are presented in Tables [TOj- 11 


4. Magnetar Bursts of unconfirmed origin 


GBM triggered on 19 bursts with magnetar-like temporal and spectral properties. Sev¬ 
eral of these bursts, however, were located in the vicinity of the Galactic Genter, which 


is a very crowded region including many conhrmed magnetars; see Table 12 for the burst 
locations. Due to the large error boxes of these bursts, we were not able to identify their 
origin unambiguously, also because no other instruments with better localization capabilities 
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recorded any of them. Figure shows the locations of all 19 bursts with their 1, 2 and 3a 
statistical error contours together with the locations of the known magnetars and magnetar 


candidates (Olausen & Kaspi 2013). All but one of the unknown bursts are in the same 


general location on the sky. The upper panel of Figure ^ displays the distribution of the 
contour centroids of the burst locations in Galactic longitude. We notice a concentration 
of 8 events between 30° and 35°. Recently, Zhou et ah (2013) reported the discovery with 


XMM-Newton of a new magnetar, 3XMM J185246.6 + 0033.7 (red diamond on Figure 21). 
The source location (/ = 33°34'45".l, 6 = —00°02'39".9) coincides with the highest density 
of unknown magnetar bursts, strongly suggesting that at least several bursts originate from 
this source (in particular those detected in 2013). 

For this sample of bursts we have also performed temporal and spectral analyses, of 
which the results are shown in Tables 12 and 13 Their analysis procedure is the same as 


we used for the bursts in Section 3G, albeit using the locations determined with the GBM 
data. 


5. Discussion Conclusions 


We have compiled here the 5-year GBM magnetar burst catalog comprising 446 events 
of which 19 are from unknown sources and 427 are distributed across seven known sources 
as shown in Table We display here the results of previously published temporal and 
spectral analyses for 357 events and include new results for the remaining 89. All these 
data are compiled here to provide a single reference catalog to facilitate large scale analyses 
of magnetar bursts for the community. The capabilities of Fermi/GEVi allow for detailed 
spectral and temporal characterizations of magnetar bursts, which can be used to further 
our understanding of these bursts. 


Figures 22 — 24 combine durations and spectral parameters for all 446 events. Each 


event is color coded per source, and the 19 unknown magnetar bursts are given the same 
color. Although SGR J1550—5418 bursts dominate the sample, we note several similarities 
in the spectral and temporal parameters across the bursts, regardless of their source. With 
respect to durations, all known source events tend to center around Tgo 100 ms (Figure!^. 
The unknown event durations tend to be shorter than most magnetar bursts, with an average 
of Tgo ~ 61 ms. However, the uncertainties on the individual durations are also considerably 
high, largely due to the low intensities of some of these bursts, hampering any dehnitive 
conclusions on the possible differences in burst durations. 


With regards to the spectral properties, the GOMP parameters of all known events are 
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very similar, with -Epeak centered at ~ 40 keV (Figure 23). The 19 unknown events have 
an arithmetic mean of F'peak 37 keV. The temperatures of the BB+BB center around 
~ 4.5 keV and ~ 15 keV (Figure [2^. The BB temperatures of the 19 unknown events have 
arithmetic means of ~ 4.3 keV and ~ 12.4 keV, very similar to the ones obtained from the 
SGR J1550—5418 bursts. For the former we hnd an average 4 ms peak flux (8—200 keV) of 
6.5 X 10“® erg/cm^/s and an average fluence of 2.3 x 10“^ erg/cm^, which places them right 
in the middle of the flux-fluence diagram of SGR J1550—5418 (see also Figure [IT|). 


Finally, the location distribution of 8 events of unknown origin strongly indicates an 
association with the recently discovered magnetar 3XMM J185246.6 + 0033.7. A dehnite 
association of these bursts is not possible, however, due to the presence of multiple magnetar 
sources in that region. 


This publication is part of the GBM/Magnetar Key Project (NASA grant NNH07ZDA001- 
GLAST, PI; G. Kouveliotou). A.G.G. was supported by an appointment to the NASA 
Postdoctoral Program at the Marshall Space Flight Genter, administered by Oak Ridge As¬ 
sociated Universities through a contract with NASA. G.K. and G.A.Y. acknowledge support 
from NASA grant NNH07ZDA001-GLAST. D.H. was supported by the Moore-Sloan Data 
Science Environment at New York University. A.L.W. acknowledges support from a Nether¬ 
lands Organization for Scientihc Research (NWO) Vidi Fellowship. A.v.K. was supported 
by the Bundesministeriums fiir Wirtschaft und Technologie (BMWi) through DLR grant 50 
OG 1101. M.v.d.K. acknowledges support from the Netherlands Organisation for Scientihc 
Research (NWO) and the Royal Netherlands Academy of Arts and Sciences (KNAW). 
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T9o (ms) T50 (ms) 


Fig. 1.— Distribution of T 90 (left panel) and T 50 (right panel) for all SGR J1550—5418 
bursts. A log-normal fit (black line) to the histograms provides /i =2.19 ± 0.02 log(ms) and 
a =0.40 ± 0.02 log(ms) for T 90 , and fi =1.71 ± 0.04 log(ms) and a =0.44 ± 0.03 log(ms) for 
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Fig. 2.— Tgo and T 50 for SGR J1550—5418 over the course of all outburst periods. The ’tail’ 
bursts in this hgure, and the following ones, are bursts detected between January 24 and 
February 24. We hnd no trends in these duration over time. 
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Fig. 3.— Evolution of the low (top panel) and high (bottom panel) BB temperatures from 
the BB+BB spectral model over time for SGR J1550—5418 bursts. For the October bursts 
we use the values for the single BB hts in both panels (red squares). There is a clear division 
between these single BB temperatures and the low temperatures of the BB+BB model. 
However, the former temperatures are similar to the hot BB temperatures of the BB+BB 
model. 
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Fig. 4.— Spectral peak as measured by the Comptonized model (top panel, i^peak) and 
OTTB model (bottom panel, kT) for SGR J1550—5418. Overall, the values are in agreement 
with each other, and they are both consistent over time. 
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Fig. 5.— Power-law index of the COMP model for SGR J1550—5418 bursts. The bursts 
from the October period are clearly different compared to the rest of the sample, conhrming 
the spectral evolution between burst active periods seen with the BB-I-BB model hts. 
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Fig. 6.— Low (left panel) and high (right panel) BB temperatures for all bursts of 
SGR J1550—5418. We £t these histograms to a Gaussian (black line) with /i = 4.55±0.05 keV 
and a = 0.71±0.05 keV for the low kT BB, and jj, = 14.96±0.16 keV and a = 2.03±0.16 keV 
for the high kT BB. 




Fig. 7.— F'peak as measured by the GOMP and OTTB models for SGR J1550—5418 bursts. 
We fit these histograms to Gaussians (black line). For the GOMP model, we find 77peak to 
have /i = 39.6 ± 0.6 keV and a =12.9 ± 0.6 keV. For the OTTB model, we find kT to have 
H = 39.8 ± 0.7 keV and a =13.7 ± 0.7 keV. 
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Fig. 8.— Power-law spectral index as measured by the COMP model for all bursts of 
SGR J1550—5418. The distribution was £t to a Gaussian (black line) with fi = —0.93 ±0.02 
and a = 0.45 ± 0.02. 
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Fig. 9.— Log N(>S) — Log (S) diagram of SGR J1550—5418 bursts. The data are fit to a 
single PL (red solid line) and a broken PL (blue dashed line). For the single PL, we only £t 
the data with S > 1 x 10“’^ erg/cm^. 
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Fig. 10.— Evolution of the 4 ms peak flux of all SGR J1550—5418 bursts as measured by 
the COMP model. The October bursts were faint compared to the other periods, but there 
is no other distinct evolution in the peak flux of the bursts over time. 
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Fig. 11.— Integrated fluence of SGR J1550—5418 bursts versus their peak flux. There is a 
strong correlation between these two parameters, and a PL £t results in a slope of 1.2 ± 0.4. 
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Fig. 12.— Tgo duration versus peak flux (top panel) and fluence (bottom panel) for all 
SGR J1550—5418 bursts. There is a marginally significant correlation between Tgo and the 
fluence, but not between Tgo and the peak flux. 
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Fig. 13.— T 50 duration versus peak flux (top panel) and fluence (bottom panel) for all 
SGR J1550—5418 bursts. There is a marginally significant correlation between T 50 and the 
peak flux, but not between T 50 and the fluence. 
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Fig. 14.— Histogram of the Tgo / T 50 ratio for SGR J1550—5418 bursts, showing that the 
distribution is broad, extending a factor of ~ 10 . 



Fig. 15.— i^peak from the COMP model versus the 4ms peak flux of SGR J1550—5418 bursts. 
We £t this correlation with a PL (dashed line) and a broken PL (solid line). 













































29 



Fig. 16.— Distribution of the phase of the peak of all the SGR J1550—5418 bursts. We do 
not hud a preferred spin phase for burst peaks, consistent with what has been found for 
other magnetars. 
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Fig. 17.— Fluence (top panel), peak flux (middle panel), and Tgo (bottom panel) of 
SGR J1550—5418 bursts compared to the spin phase. The gray dots represent each indi¬ 
vidual burst, and the black dots represent the weighted average of these gray dots grouped 
by phase. The size of the black dots represent how many bursts are in each group. The blue 
diamonds indicate the 0.2—10 keV pulse prohle (Lin et ah 2012). There are no correlations 
between the fluence, peak flux or Tgo and the spin phase. 
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Fig. 18.— Folded lightcurve of all bursts from SGR J1550—5418 in 3 channels: 8—20 keV 
(blue), 20—50 keV (red), and 50—200 keV (yellow). The sum of all three channels is shown in 
green. There is a possible peak at the bin that spans phase 0.15—0.20. This peak, however, 
is not signihcant (see Section 3.1.5). 
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Fig. 19.— Folded lightcurve of all bursts from SGR J1550—5418 without the ten brightest 
ones (top panel), and folded light curve of only the ten brightest ones; color coding the same 


as in Figure 18 The peak at phase bin 0.15—0.20 in Figure 18 is not apparent in the top 


panel, but clearly present in the bottom panel. 
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Fig. 20.— Log N(>S) — Log (S) diagram of SGR J0501+4516. The data are fit to a single 
PL (red solid line) and a broken PL (blue dashed line). For the single PL, we only £t the 
data with S > 1 x 10“^ erg/cm^. 
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Galactic Longitude (deg) 


Fig. 21.— Locations in galactic coordinates of magnetar-like bnrsts with nnconfirmed origin, 
displayed as 1, 2, and 3a statistical error contonrs (bine lines). The known magnetars and 
magnetar candidates are shown as black diamonds, and the magnetar 3XMM J185246.6 + 
0033.7 (Zhon et al. 2013) is highlighted in red. The histogram in the top panel is for the 


centroids of the GBM location contours. 
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Fig. 22.— Distributions of the T 90 and T 50 durations for all magnetar bursts in this catalog. 
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Fig. 23.— Distributions of F^peak and the spectral index of the COMP model for all magnetar 
bursts in this catalog. 
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Fig. 24.— Distributions of the low- and high-temperature BB for all magnetar bursts in this 
catalog. 



















































- 38 - 


Table 1. Summary of GBM magnetar bursts 


Source 

Burst active periods 

Number of bursts 

with TTE data 

SGRJ1550-5418 

Oct 2008 - Apr 2009 

386 

SGRJ0501+4516 

Aug/Sep 2008 

29 

IE 1841-045 

Feb - Jul 2011 

6 

SGRJ0418+5729 

Jun 2009 

2 

SGR1806-20 

Mar 2010 

1 

SGRJ1822.3-1606 

Jul 2011 

1 

AXP4U 0142+61 

Jul 2011 

1 

AXP IE 2259+586 

Aug 2011 

1 

unknown 
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Table 2. Durations, phases, peak fluxes, and fluences for SGR J1550—5418 bursts 
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Table 6. Durations, peak fluxes, and fluences of SGR J0501+4516 bursts 
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Table 7. Time-integrated spectral modeling results for SGR J0501-1-4516 bursts 
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Table 9. Time-integrated spectral modeling results for IE 1841—045 bursts 
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